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Abstract 

"7 + Jet" events, based on the qq ^ g + J and qg ~> q + J subprocesses, are proposed for jet energy 
scale setting and hadron calorimeter calibration at LHC energies. General features and selection criteria of 
"7 + Jet" events that would provide a good Pf^ — Pj-^*^* balance are described. CMS detector geometry is 
taken as the basement. 



1. INTRODUCTION 



Setting an absolute scale for a hadronic calorimeter (HCAL) is an important task for many 
of pp and pp collider experiments (see e.g. [1, 2]). There is a number of ATLAS and 
CMS publications on this subject where the application of different physical processes (Uke 
"Z°/7 + Jet" and others) is discussed ([3-10]). 

This paper is the first part of a set of publications on a detailed study of this problem. 

The main goal of this work is to find out the selection criteria for "7 + Jet" events that 
would lead to the most precise determination of transverse momentum of a jet, i.e. Pj'^'^*, via 
assigning a photon to a jet. We shall present here the results of event generation by using 
PYTHIA 5.7 [11]. Further development, based on the simulation of detector response with 
GEANT detector simulation packages will be presented in one of our next papers. 

It should be also noted that here we consider the case of a low luminosity L = 
j^q33 jj^j^j jj, gjjn q^jl■g sufficient to use much more restrictive cuts as well as new 

physical variables and, correspondingly, cuts on them (cluster suppression criterion, for in- 
stance) in comparison with those used in the previous experiments and to obtain a set of more 
clean "7 + Jet" events. 

The Section 2 is an introduction to the problem. General features of the "j+Jet" pro- 
cesses at LHC energies, that will be explored in this article and later papers, are presented 
here. In Section 2 we review possible sources of the Pt"^ and p^-^®* disbalance in the final 
state and the ways of selecting those events where this disbalance has a minimal value. 

In Section 3.1 we present definitions of components of different objects that enter 
the balance equation illustrating the conservation law of the total Pt in any event. 

Section 3.2 describes the criteria we have chosen to select "7 + Jet" events for the 
calibration procedure. "Cluster" (or mini-jet) suppression criterion iPi^Q'^f^) which has not 
been used in previous experiments is introduced here. Its important role will be illustrated in 
the following papers [12-15]. These clusters have a physical sense as a part of another new 
experimentally measurable quantity introduced here for the first time, namely, the sum of the 
Pt of all particles detectable in the [77 1 < 5 region which are out of the "7 + Jet" system 
(denoted as Pt""*). 

Another new thing here is an introduction of a new physical object, named as an 
"isolated jet", i.e. the jet that does not have any noticeable Pt activity in some ring in — ^ 
space around it. In other words we will select some class of events having a total Pt activity 
inside the ring, AP = 0.3, around this "isolated jet" within 2 — 5% of jet Pt- In the following 
paper it will be shown that the number of events with such a clean topological structure would 
not be small at LHC energies (mainly due to the growth of luminosity). 

Since the calibration is rather a practical than an academical task in all the following 
Sections, we present the values of rates for strict and weak cut values because their choice 
would be a matter of step-by-step collected statistics. 

The justification of the variables and cuts introduced in Section 3 can be found in our 

papers [13-15]. In [15] we present the estimation of the efficiency of background suppres- 
sion, that is, finally, the main guideline to establish the selection rules. 

Section 4 will be devoted to the estimation of non-detectable neutrino contribution to 
P('^^* as well as to studying the influence of the I77I > 5 region not covered by calorimeters or 
other detectors (that is the main source of Pj™'"* =-^t) on the total Pt balance in the event. 
The correlation of the upper cut on p^"**"* with a mean value of Pt of neutrinos belonging to 
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the jet P(, i.e. {Ptf^^}, will be considered here. 

Since the results presented here have been obtained with PYTHIA simulation, we 
are planning to carry out the analogous estimation in the next papers but with another event 
generator. 

2. GENERALITIES OF "7 + Jet" PROCESSES 
2.1 Leading order picture 

The idea of a hadronic calorimeter (HCAL) calibration by physical process "pp ^ 7 + Jet" 
was realized many times in different experiments (see recent papers [1,2] and refs. there). 
It is based on the parton picture where two partons (qq or qg), supposed to be moving in 
different colliding nucleons with zero transverse momenta (with respect to the beam line), 
produce a photon, called a direct one, and a parton with balanced transverse momentum 

-* part 7 

Pt = —Pt ■ This picture corresponds to the leading order (LO) Feynman diagrams 
shown in Fig. 1 for the "Compton-like" process 

+ 1 (la) 

and for "annihilation" process 

qq-> g + J, (16) 

respectively. The Pt of the "7+parton" system produced in the final state should be equal to 
zero, i.e. 

-* 7+r)art ~* part 

Pt =Pt +Pt =0. (2) 

So, in this case one could expect that with a reasonable precision the transverse momentum of 
the jet produced by the final state parton (q or g) will be close in magnitude to the transverse 

-* Jet -* 7 

momentum of the final state photon, i.e. Pt ~ —Pt ■ 




Fig. 1 : Some of the leading order Feynman diagrams for direct photon production. 

It allows one to carry out the calibration of HCAL in the experiments with a well 
calibrated electromagnetic calorimeter (ECAL). To put it simpler, to a part of jet transverse 
energy E^^* deposited in HCAL we can assign the value of the difference between the value 
of the transverse energy deposited in ECAL in the photon direction, i.e. E^, and the value of 
the transverse energy deposited in ECAL in the jet direction. 

2.2 Initial state radiation. 

Since we believe in the perturbation theory, the leading order (LO) picture, described above, 
is expected to be dominant in determination of the main contribution to the cross section. 
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The Next-to-Leading Order (NLO) approximation (see some of the NLO diagrams in Fig. 2) 
introduces some deviations from a rather straightforward LO-motivated idea of caHbration. 
Thus, as it is seen from Fig. 2, a gluon radiated in the initial state (ISR) can have its own 
non-zero transverse momentum p^^^^ ^ Q. It leads to the non-zero transverse momenta of 
partons that appear in the initial state of fundamental 2^2 QCD subprocesses (la) and 
(lb). As a result of the transverse momentum conservation, a disbalance of the transverse 
momentum of a photon P^^ and of a parton p^P'*' ' produced in the fundamental 2^2 pro- 
cess 5 + 6^7-1-8, shown in Figs. 2 and 3 (and thus, finally, of a jet produced by this parton), 
wiU take place. 





Fig. 2: Some of Feynman diagrams of direct photon production including gluon radiation in the initial state. 



^6 



We shall choose the modules of a vector sum of transverse momentum vectors Pt and 
Pt of incoming (into 2^2 fundamental QCD subprocesses 5 -|- 6 — > 7 -|- 8) partons, shown 
on lines 5 and 6 in Fig. 2, as well as the sum of their modules as two quantitative measures 



o5+6 



t =\p; + pi\, Pt56 = IP* V iPt*^ 



to estimate the Pt disbalance caused by ISR. The modules of the vector sum 



\Pt +Pt 



(4) 



will be used as an estimator of the final state Pt disbalance in the "7 -|- Jet" system. 

The numeration notations in these Feynman diagrams as well as in formulae (3) and 
(4) are chosen to be in correspondence with those used in PYTHIA for describing parton- 
parton subprocess, displayed schematically in Fig. 3. The "ISR" block describes the initial 
state radiation process that can take place before the fundamental hard 2 — > 2 process. 
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Fig. 3: PYTHIA "diagram" of 2 2 process (5+6-+7+8) with a block (3+4-»5+6) of initial state radiation (ISR). 



2.3 Primordial parton kr effect. 

A possible non-zero value of the intrinsic parton velocity inside a colliding proton may be 

another source of the Pt'^ and p^P"-^* disbalance in the final state . Its reasonable value is 
supposed to lead to the value of kx < 1-0 GeV/c. It should be noted that sometimes in the 
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literature the summarized effect of ISR and of the intrinsic parton transverse momentum is 
denoted by kx- Here we follow the approach used in PYTHIA where these two sources of 
disbalance are treated separately and switched on by different keys (MSTP(61) for ISR and 
PARP(91), PARP(93) and MSTP(91) for kr). Below we shall keep the value of kr to be 
fixed by PYTHIA default value (fcji)=0.44 GeV/c. Its possible variation influence on Pi'' 
and Pj-^*^* disbalance will be discussed in detail in our following paper [15]. The general 
conclusion from there is that any variation of kr within reasonable boundaries (as well as 
beyond them) does not produce a large effect in the case when the initial state radiation is 
switched on. The last one gives a dominant contribution. 

3. CHOICE OF MEASURABLE PHYSICAL VARIABLES FOR "7 + Jet" PRO- 
CESS AND OF CUTS FOR BACKGROUND REDUCTION. 

Another than (la) and (lb) QCD processes with large cross sections, being by orders of 
magnitude higher than the "7 + Jet" cross section, can also contain high Pt photons and 
jets in final states. So, we face the problem of signal "7 + Jet" events selection out from the 
large QCD background. Here we shall discuss the choice of physical variables that would be 
useful under some cuts on their values to select the desirable processes with direct photons 
("7*''") production from the background events. The possible 7*'' candidate may originate 
from the TT—, t]—, oj—, if — meson decays or, may be, from a bremsstrahlung photon. 

We suppose the ECAL size to be limited by I77I < 2.61 and HCAL to consist of the 
barrel (HB), end-cap (HE) and forward (HP) parts and to be limited by \ri\ < 5.0 (CMS ge- 
ometry), where rj = —0.5 In {tan {0/2)) is a pseudorapidity defined through a polar angle 9 
counted from the beam line. In the plane transverse to the beam Une the azimuthal angle (f) 

Jet ~* 7 

defines the directions of Pt and Pt . 



3.1 Introduction of some new measurable physical observables and Pt balance equa- 
tion. 

In — > 7 + Jet + X events, we are going to study, the main physical object will be a high 
Pt jet, to be detected in the 1 77] < 5 region, and a direct photon, registered by ECAL up to 
< 2.61. In these events there will be a set of particles mainly caused by beam remnants, 
i.e. by spectator partons fragmentation, that are flying in the direction of a non-instrumented 
forward part {\rj\ > 5) of the detector. Let us denote the total transverse momentum of these 
non-observable particles as 

E P:^!"'""- (5) 

ie|r)|>5 

Among the particles with < 5 there could be also neutrinos. Their total momentum 
will be denoted as 

E Pt\.)^Pti.y (6) 
ie|??|<5 

The sum of transverse momenta of these two kinds of non-detectable particles will be denoted 
asPt™"": 

miss -* r7>5 

Pt =Pt^,)+Pt^ . (7) 
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A high energetic jet can also contain neutrinos that may carry some part of the total 
jet energy needed to be estimated from simulation. From the total jet transverse momentum 

Jet 

Pt we shall separate the part that can be measured in the detector, i.e. in the ECAL+HCAL 

jet 

and muon systems. Let us denote this part as Pt (small "j"!). So, we shall present the total 

^ Jet 

jet transverse momentum Pt as a sum of three parts: 

^ Jet 

1. Pt{i,) ■, containing the contribution of neutrinos that belong to the jet, i.e., a non- 
detectable part of jet Pt- 

-> Jet r — ^ i 
Jet -* Jet 

2. Pt(i_t) J containing the contribution of jet muons to Pt . These muons can give a 
weak signal in the calorimeters but their energy can be measured in the muon system (in the 
region of |r/| < 2.4 in the case of CMS geometry): 

-» Jet ^ — ^ i 
iGJet 

Let us mention that due to the absence of the muon system and tracker beyond the 
77I < 2.4 region, there exists a part of p^'^"^* caused by muons with \ri\ > 2.4 denoted as 
Pt'llf\ri\>2 4)- This part can be considered in some sense as an analog of Ptf^) since the only 
trace of its presence would be weak MIP signals in ECAL and HCAL. 

As for both points 1 and 2, let us say in advance that the estimation of the average 
values of the neutrino and muon contributions to Pt '^'^* (see Section 4 and also Tables 1-8 
of Appendix) has shown that they are quite negligible: about 0.35% of {Pt'^'^*')aU is due to 
neutrinos and about 0.25% of {Pt'^^'^)aii — to muons, where all means averaging over all 
events including those without neutrinos in jets. 

jet Jet 

3. And finally, by means of Pt we denote the part of Pt which includes all 
detectable particles of the jet , i.e.: 

^ 7ci * Jet * Jet 

Pt =PtiHCAL+ECAL) + Ptif.), Iv^l < 2A (10) 

Thus, we can write (for the general case of 77 values) 

Jet jet Jet -» Jet 

Pt = Pt + Pt(iy) + Pt(^,|,,f |>2.4) (11) 

In the case of pp ^ j + Jet + X events the particles detected in the jTyj < 5 region 
can originate from the fundamental subprocesses (la) and (lb), that may be caused by LO 
diagrams, shown in Fig. 2 as well as by NLO diagrams (like those in Fig. 3 that include ISR) 
and also from the "underlying" event, of course. 

As it was already mentioned in Section 2, the final states of the fundamental subpro- 
cesses (la) and (lb) may contain additional jets due to the ISR and final state radiation (FSR) 
caused by the higher order QCD corrections to Feynman diagrams given in Fig. 1. To real- 
ize the calibration idea (see Section 2.1), we need "in situ" selection of events with a good 

-•7 jet 

balance of Pt and the Pt part, measurable in the detector It means that to make a reason- 
able simulation, we need to have a selected set of events with a small Pt(u) (and, thus, small 

' We shall consider the issue of small Pt charged pailicles contribution into the total jet Pt while discussing the 
results of the full GEANT simulation (with account of the magnetic field effect) in our forthcoming papers. 
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p^miss-^ as a model and we also have to find a way to select events without additional jets or 
with jets suppressed to the level of very small Pt mini-jets or clusters. 

For any event let us separate the particles in the < 5 region into a "7 + Jet" sys- 
tem (here "Jet" denotes the jet with the highest Pi > 30 GeV/c) having the total transverse 

^ 7+ Jet 

momentum P^ (see (4)) that may be different from: 

Pt =Pt +Pt , (12) 

in the case of neutrino presence in a jet, and a system of aU other (O) particles in the |?7| < 5 
region beyond the "7 -|- Jet" system. The total transverse momentum of this system will be 
denoted as P^*^ and it is a sum of Pt of additional mini-jets (or clusters) as well as of Pf of 
single hadrons, photons and leptons in the |??| < 5 region. Since neutrinos are present among 

-. -. Jet 

these leptons, then the difference of Pt{v) and Pt(v) gives us the value of the transverse mo- 
mentum 

PtZ)=Pti.)-Ptl^, Kl<5 (13) 

carried out by neutrinos which do not belong to the jet but are contained in the |7/| < 5 region. 
Let us denote a part of Pt , that can, in principle, be measured in the detector, by 

-t out -* out I ^ 

Pt ■ Thus, Pt is a sum of Pt of other mini-jets or clusters (with Pt smaller than 
Pt "'''*) and of Pt of single hadrons (h), photons (7) and electrons (e) with < 5 and muons 
(11) with |?7^| < 2.4. Below for simphcity these mini-jets and clusters wiU be called just 

out 

"clusters". So, Pt is the following sum (h, 7, e, ^ Jet): 

-* out dust sing nondir -* -* O 

Pt = Pt + Pt(h) + Pt{-/) + Pt(e) + Pt{ix,\rit'\<2A)l 1^1 < 5 (14) 

And thus, finally, we have: 

^out ->0 

Pt = Pt + Pt(u) + Pt(f_i^\jjii\>2A)- (15) 

With these notations the conservation law for the Pt component of the whole "7 + Jet" 
(where 7 is a direct photon) event is: 

^7 Jet ^))>5 

Pt +Pt +Pt +Pt =0 (16) 
with last three terms defined correspondingly by (11), (15) and (5). 

3.2 Definition of selection cuts for physical variables. 

1. We select the events with one jet and one photon with 

Pt'' > 40 GeV/c ; Pt'"'' > 30 GeV/c. (17) 

For most of our applications the jet is defined according to PYTHIA jetfinding algorithm 
LUCELL. The jet cone radius R in — space , counted from the jet initiator cell (ic), is 

taken as Ric = ((A?])^ + (A(/))^)-^/^ = 0.7 Below we shall also use the value of the jet 
radius, counted from the center of gravity (gc) of the jet, i.e. Rgc. The comparison with UAl 
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jetfinding algorithm (taken from CMSJET program of fast simulation [16]) is presented in 
[13, 14]. 

2. To suppress the background processes, only the events with "isolated" photons are taken. 
To do this, we restrict: 

a) the value of the scalar sum of Pt of hadrons and other particles surrounding a pho- 
ton within a cone of Rj^^i = {{Ar])'^ + {A^)'^y/'^ = 0.7 ("absolute isolation cut") 

J2Pt' = Pt''"' <PtcuT; (18) 

b) the value of a fraction ("relative isolation cut") 

Y,Pt'/Pt''^e'^<eluT-, (19) 
ieR 

c) we accept only the events having no charged tracks (particles) with Pt > 1 GeV/ c 
within i?7soi ^^^^ around the photon candidate. 

3. To be consistent with the application condition of the NLO formulae, one should avoid an 
infrared dangerous region and take care of Pt population in the region close to a photon (see 
[21-23]). In accordance with [22] we also restrict the scalar sum of Pt of particles around a 
photon within a cone of a smaller radius Rsingi = 0.175 = 

Due to this cut, 

= p^smsi < 2 GeF/c, {i^'y-dir). (20) 

an "isolated" photon with high Pt also becomes a "single" one within an area of 8 towers (of 
0.087x0.087 size according to CMS geometry) which surround the tower bitted by it (analog 
of 3 X 3 tower window algorithm). 

4. We also consider the structure of every event with the photon candidate at a more pre- 
cise level of 5x5 crystal cells window (size of one CMS HCAL tower) with a cell size of 

0. 0175x0.0175. To suppress the background events with photons resulting from high ener- 
getic 7r°— , T]—, uj~ and Kg— mesons, we apply in addition the following cut: 

a) the ECAL signal can be considered as a candidate to be a direct photon if it fits 
inside the 3x3 ECAL crystal cell window with the highest Pt 7/e in the center; 

b) the value of a scalar sum of Pt (p^ of stable particles in the 5x5 crystal cell 
window in the region out of a smaller 3x3 crystal cell window (typical size of photon shower 
in ECAL found from GEANT simulation with CMSIM package), having the cell with the di- 
rect photon candidate (the largest 7/e) as the central one, should be restricted by 1 GeV/ c, 

1. e. 

p^sum < ^ GeV/c; (21) 

c) we require the absence of a high Pj hadron in this 5x5 crystal cell window (that 
means an imposing of an upper cut on the HCAL signal at least in the one-tower area) around 
the direct photon: 

p^hadr ^ 5 GeV/c. (22) 

We can not reduce this value to, for example, 2-3 GeV/c, because a hadron with Pt 
below 2-3 GeV/c deposits most of its energy in ECAL and may not reveal itself in HCAL. 
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^ Jet ^ 7 

5. We select the events with the vector Pt being "back-to-back" to the vector Pt (in the 
plane transverse to the beam Une) within A(j!) defined by equation: 

0(^,,et) = 180° ± A,^ (A,/. = 15°, 10°, 5°) (23) 
(5° is a size of one CMS HCAL tower in cj)) for the following definition of the angle 0(7,jet) 

P,^P/^* = P,^P/«* . COs(.^(^,,et)), With PjT = |P/|, P/^* = |P/'*|. 

6. The initial state radiation (ISR) manifests itself as some final state cluster or mini-jet ac- 
tivity. To suppress it, we impose a new cut condition that was not used earlier in previous 
experiments: we choose the events that do not have any other jet- like or cluster high Pt activ- 
ity, i.e. Pj°'"** (taking the cluster cone Rciust {v^ 0) = 0-7), being higher than some threshold 
PtcuT value, i.e. we select the events with 

PI dust ^ D clust f^A\ 
t S i^tcUT^ U"*) 

where clusters are found by one and the same jetfinder used to find the jet in the event. 

7. We limit the value of modulus of the vector sum of Pt of all particles, except those in the 
"7 + Jet" system, that fit into the region covered by ECAL and HCAL (i.e. the cells "beyond 
the jet and photon" regions): 



E 



i^Jet^^—dir 



Pt"""' < PtcliT. 1^1 < 5 (25) 



The importance of PtcuT PtciiT parameters to reduce the background will be demon- 
strated in the forthcoming papers [13-15]. 

Below the selection cuts 1-7 will be referred as "Selection 1". The last two of them, 
6 and 7, are new criteria, not used in previous experiments. In addition to them one more new 
object, named an "isolated jet", will be introduced. 

8. To do this we also involve a new requirement of "jet isolation", i.e. the presence of a 
"clean enough" (in the sense of limited Pt activity) region inside the ring (of AR = 0.3 size) 
around the jet. Following this picture we restrict the value of the ratio of the scalar sum of 
transverse momenta of particles belonging to this ring, i.e. 

= e^"* < 2%, where Pt""^ = ^ \Pt"\. (26) 

ieO-7<R<l 

The set of events that pass under the cuts 1-8 will be called as "Selection 2". 

9. In the following "Selection 3" we shall keep only those events in which one and the 
same jet (i.e. up to good accuracy having the same values of Pt^^\ W^* and A0) is found 
simultaneously by every of two jetfinders used here: UAl and LUCELL. For these jets (and 
also clusters) we require the following conditions: 

P/'^* > 30 GeV/c, Pt < PtcuT, < 15°(10°, 5°), e-''^* < 2% (27) 

The exact values of cut parameters, i.e. PtcuT> ^cut' Ptcur^ PtcuT' will be 
specified below, since they may be different, for instance, for various Pt^-intervals (being 
more loose for higher Pt''). 
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10. As we have already mentioned in Section 3.1, one can expect reasonable results of the 
calibration procedure modeling only by using a set of selected events with a small value of 
p^miss ^jgQ j^gg following cut: 

P> miss ^ r> miss 
t S "tCUT- 

Due to this reason in Section 4, we shall study the influence of parameter on the 

selection of events with a reduced value of Ptfu)- The aim of the event selection with a 
small value of Ptf^^ is quite obvious: we need a set of events with a reduced value of Pj-^"^* 
uncertainty due to possible presence of a non-detectable neutrino contribution to a jet. 

To conclude this section, let us rewrite the basic Pf -balance equation (16) of the previ- 
ous section by means of notations introduced here in the form more suitable to present results 
in further papers [12-15]. For this purpose we shall rewrite equation (16) in the following 
scalar form: 

pi _ p -fe* 

* = (l-cosA(/.) + Pt(0 + 77> 5)/Pt\ (29) 

where Pt{0 + rj > h) = {Pt° + P*''''^^^) • n-^"^ with iV"^ = p/''*/P/''\ 

As it will be shown in [12-14], the first term in the right-hand part of the equation 
(29) is negligibly small and tends to decrease more with a growth of the energy. So, the main 
source of the Pt disbalance in the "7 + Jef system is a term Pt(0 + 77 > 5) /Pt^. 



4. ESTIMATION OF NON-DETECTABLE PART OF Pt ' 

In Section 3.1 we have separated the transverse momentum of the jet, i.e. Pt'^'^^ , into two 
parts: a detectable one Pt-*"^* and a non-measurable part, consisting of Ptf^^ (see (8)) and 
Ptflf\ri\>2 4) (s^^ (11)) 0- In the same way we have done analogous separation according 
to equation (15) of the transverse momentum of other particles, i.e. Pt*^, excluding direct 
photon (or candidate to be detected as a direct photon), into detectable part p^""* and non- 
measurable part consisting of Pt(^) (see (13)) and Pt(^j,,|>2.4) (see (15)). 

Here we present an estimation of averaged values of transverse momenta of the total 
Pt'''^* carried out by non-detectable particles. For this aim we use a bank of the signal "7 + 
Jef events generated for three intervals of Pt'' with the restrictions (17) - (24) and the 
following cuts are fixed as follows: 

Pt}^°Jj, ^20 GeV/c, eluj. = 15%, A0 = 15°, PtcuT = 30 GeV/c. (30) 

No restriction for the Pj°"* value was done. The results of analysis of these events are pre- 
sented in Figs. 4 and 5. 

The first row of Fig. ^ contains p^""*^ spectra in the "7 + Jef events for different 
Pt^ intervals. Their practical independence (up to the good accuracy) on Pt^ is clearly seen. 

In the second row of Fig. Q we present the spectra of p^""''* for the events (denoted 
as Ptf^^ > 0) having a non-zero Ptf^^ component in P*'^''*. For these figures the Pj™"**^ 
spectrum dependence on the direct photon Pt'' (that is equal, approximately, to Pt"^*^*) is seen. 

-Firstly we shall consider the case of switched off decays of and mesons (according to the PYTHIA 
default agreement, vr^ and mesons are stable). 
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So, the spectra tails as well as the mean values are shifting to a large region with p^-''^* 

growth. At the same time a peak position remains in the region of p^™*^'* < 5 GeV/c. From 
the comparison of the number of entries in the second row plots of Fig. ^ with those in the 
first row it can be concluded that the part of events with the jet having the non-zero neutrinos 
contribution {Ptf^^ > 0) has the same size of about 3.3% in all Pt'' intervals. 

The same spectra of p^™"*'' for events with Ptf^^ > show what amount of these 
events would remain after imposing a cut on p^""^^ in every Pt'' interval. The important 
thing here is that the reduction of the number of events with Ptf^* > in every Pf ^ interval 

leads to reducing the mean value of Ptf^*, i.e. a value averaged over all collected events. 
This value, found from PYTHIA generation, serves as a model correction and it has to 
be estimated for the proper determination of the total Pj'^'^* from the measurable part P*-''^*: 

p^Jet ^ p^jet ^ ^ {Ptl:^)aU events- 

The effect of general PtcuT imposing in each event of our sample is shown in the 
third row of Fig. ^ The upper cut PtcuT — 1000 GeV/c, as it is seen from the second row 
pictures, means the absence of any upper limit on Ptf^y The most important information 
that the value of the neutrinos Pt inside the jet, being averaged over all events, can reach 
the value of Ptf^;) w 1 GeV/c at PJ > 300 GeV/c comes from the right-hand plot of 
the third row in Fig. 4. From the comparison of the plots from the second row with the 
corresponding plots ^ from the third row we see that the first cut Ptcur — 20 GeV/c for the 
first 40 < Pt^ < 50 GeV/c interval reduces the number of entries by less than 0.5% and the 
mean value of Ptjv* — by less than 10%. A more restrictive cut PtcuT — 5 GeV/c reduces 
the value of {Ptf^*) by three times and leads to approximate twofold drop of the number of 
events. 

From these Figures we see that for the interval 300 < Pt'' < 360 the number of events 
with jets including neutrino (second row) is about 3.3% (Entries=3001) of the total number 
of the generated "7 + Jet" events (Entries=89986). A very restrictive P(™^'*|,=5 GeV/c cut 
leads to the reduction factor for (Ptf^^) of about 50. As it is seen from the plot in the bottom 
right-hand corner of Fig. Q the Ptf^) spectrum for the remaining events (Entries=57475) 
finishes at Pt'^^ = 10 GeV/c and sharply peaks at Pt"^^* = 0. The averaged value of Pt'^^^ 
under this peak is equal to 0.022 GeV/c. So, with this cut on p^™*"* the neutrinos give a 
negligible contribution to Pt"''^*. 

At the same time we see that application of the moderate cut Pt^cuT — 10 GeV/c 
for 300 < Pt'^'^* < 360 GeV/c interval strongly reduces (by 20 times) the mean value of 
Ptf^) (from 1 GeV/c down to < Ptf^) >= 0.05 GeV/c) at about 10% reduction of the total 
number of events in this p^-^'^* (or Pt'') interval. 

In the case of 100 < p^'^'^* < 120 GeV/c interval, as we see from the third row of 
Fig. the same cut Ptcur 

= 10 GeV^/c reduces the mean value of Ptfy^ by 5 times (from 

0.5 GeV/c down to < Pt(^^ >= 0.09 GeV/c) with the same 10% reduction of the total 
number of events. 

It should be noted that in the less dangerous (from the point of view of the size of 



^That includes the values of Pt^lf^ ^nd the coiTesponding number of entries remained after Pt^}f^ imposing 
as well as the mean value of Pt'^^y denoted as "Mean" (being equal to an averaged {Pt'^^^) value over the number 
of the remained entries. 



10 



40 <Pt^< 50 GeV/c 



100 <Pt^< 120 GeV/c 



300 <Pt^< 360 GeV/c 






Fig. 4: a) spectra in all events; b) Pj""'"* spectra in events having jets with non-zero Pt neutrinos, i.e. 

Pt/fj* > 0; c) -Pt/^N* spectra behavior for different values of PtcuT values in various Pt^''*{Ki Pt'^) intervals. 
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40<Pt^<50 GeV/c 



100<Pf< 120 GeV/c 



300<Pt^<360 GeV/c 





Fig. 5: a) p^™*"' spectra in all events; b) p^™*"^ spectra in events having jets with non-zero Pt neutrinos, i.e. 
Pt/jj* > 0; c) ft^j)* spectra behavior for different values of PtcuT values in various Pt-^'=*(RJ Pt^) intervals. 
Ji'^'^— decays are allowed inside the solenoid of ii = 129 cm and L = 317 cm. 
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neutrino Pt content in a jet) 40 < Pj-''^* < 50 GeV/c interval we have already a very small 
mean value of Ptf^^ equal to 0.152 GeV/c without imposing any PtcuT- 

The analogous (to neutrino) situation holds for Ptfl^^ contribution. 

The detailed information about the values of non-detectable Ptfu)^ averaged over all 
events (no cut on p^™-^^ was used), as well as about mean values of muons Pt from jet 
Ptf^), is presented in Tables 1-8 of Appendix for the sample of events with jets which are 
completely contained in the barrel region of HCAL (\r]^'^^\ < 1.4, "HB-events", see [12]). 
In these tables the ratio of the number of events with non-zero Ptf^^ to the total number of 
events is denoted by R^^^nt ^aAo of the number of events with non-zero Pt j'^* to the 

total number of events is denoted by R'^^/nt ■ The Pt quantity in events with Pt f^j > is 
denoted in these tables as Pt"''iyet ™d is given there for four Pt'* intervals and other Pt'cijT 
values than in the second row plots of Figs. 4 and 5. From Tables 1-8 we see that the averaged 
of value being calculated by using of only the events with Ptf^) > 0, i.e. {PtZ^Jet)' 

is about 7 GeV/c for 40 < Pt'^ < 50 GeV/c interval. It increases to about 32 GeV/c for 
300 < Pt'^ < 360 GeV/c. It should be noted that the averaged values of the modulus of Ptf^^ 
(see formula (8)) presented in the second and third lines of Tables 1-6 from the Appendix, 
coincide with the averaged values of up to the three digits due to practical collinearity of 

Jet jet 

Pt and Pt vectors, i.e. <Pt(^) >=<A^>. 

Tables 1-8 contain the Unes with an additional information on the numbers of the 
"7 + Jei" events with the jets produced by c— and b— quarks, i.e. N event (c) and Nevent{p), 
given for the integrated luminosity Lint — 3 fh^"^, respectively. There are also lines that 
show a ratio ("29sub/air) of the number of the events caused by gluonic ("Compton-like") 
subprocess (la) to the number of events due to the sum of (la) and (lb) subprocesses. 

Below follows the line containing the averaged values of the jet radius <Rjet>- 
The value of the difference of the jet measurable part transverse momentum Pt-**^* and 
of the total jet Pt^^*, averaged over all events, i.e. {Pt'''^^ — Pt"*^^), is presented in second 
lines of Tables 1-8 of Appendix in GeV/c units. This value has a sense of correction 
that should be applied to Pf-''^* in order to take into account Pt carried off by non-detectable 
particles. 

It was already mentioned in Introduction that we are planning to carry out a more 
detailed analysis basing on GEANT package and taking into account weak decays of vr^ and 
mesons. To have an idea what changes can be expected, we shall consider only the case 
of allowed decays (as the main source of neutrinos and muons). The averaged values 
of Pt'll^^ for different Pt^-intervals with switched on decays are given in Fig. ^ with the 
same meaning of all notations. Here decays are allowed inside the solenoid volume with 
the barrel radius Rb = 129 cm and the distance from the interaction vertex to End-cap along 
Z-axis L — 317 cm (CMS geometry). 

From this Figure we see that in the case of allowed decays, the Pj"""'* spectrum 
for all events (compare the first rows in Fig. 4 and Fig. 5) practically does not change with 
p/'=*(« Pt'^) as well as the mean value of Pf"""", At the same time the p^""*** spectra 
change for events that contain neutrinos in the jet (second row of Fig. 5) quite noticeably. It 
should be noted that the amount of such events grows up to 10% as compared with 3% in 
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the case considered in Fig. 4, but the mean values of p^""-^" do not grow so strongly with 
Pt'' as it was seen in Fig. 4. The mean value of p^"^^^" changes only twice from 6.0 GeV/c 
for the interval 40 < Pt'^ < 50 to 11.4 GeV/c for the interval 300 < P{^ < 360. Now 
we compare the third row pictures in Figs. 4 and 5. We see that the mean value of Ptf^^, 
carried out by neutrinos of the jet grows up from the value of about {Ptfv)) ~ 0-^ GeV/cdX 
40 < < 50 GeV/c to the value of {Ptf^^) « 2.2 GeV/c at 300 < P{^ < 360 GeV/c, i.e. 
it changes as 2% 0.7% of Pi From the same pictures of Fig. 5 we see that the general 
cut PtcuT = 20 GeV/c would reduce the contribution of neutrinos into Pt"^^* to the value 
(Ptfj)*) » 0.5 GeV/c in all P/^* intervals, while the cut Pt^^'^ = 10 GeV/c would lead to 
the value of {Pif^)) « 0.20 - 0.26 GeV/c. (that is quite acceptable) with k. 9% reduction 
of the event number. 

5. SUMMARY 

The possibility of the jet energy scale setting and hadron calorimeter calibration at LHC 
energies by using the "7 + Jet" process is considered for the case of low luminosity 

(1033 cm^s-i). 

The initial state radiation (ISR) as the main source of the Pt disbalance in the "7 + 

Jet" system is discussed. New variables that enter the P( -balance equation are considered. 
The new cuts (see Section 3) Ptcur ^nd PtcuT (i^^ addition to the cut on (/^(■yjet) angle 
used previously in other experiments) as well as a new object of "isolated jet" are introduced 
here. The consequences of their variation and the choice of their most preferable values to 
select the events with a good Pt'^ and P*-''^* balance will be discussed in [12-15]. 

The values of the non-detectable part of Pj -''^* caused by neutrinos are estimated for 
different p^"^^* intervals. It is found that Pt'' and p^'^'^* balance can be influenced by neutrino 
energy leakage from the jet. The PtcuT = 10 GeV/c is proved to be sufficient for reducing 
the jet energy leakage caused by neutrinos to an acceptable level with about of 9% loss of 
events. 

After the detailed study of neutrinos and muons contribution to Pf'''^* done here, the 
following papers will be concentrated on the contribution of hadrons, photons and electrons. 

We would like to emphasize once more that the values of selection cuts, given here, 
are not a dogma for us. Our aim is to present an estimation of a number of events that can 
be selected in some unit of time chosen here as one month of LHC continuous operation (i.e. 
3000 = 3 fb~^). In future calibration "in situ" one can collect these "7 -\- Jet" events 
and in parallel classify them according to different Selection criteria (1,2 and 3; see Section 
3.2) for increasing the degree of accuracy. 
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APPENDIX 

40 < Pt^ < 50 GeV/c 



Table 1: Selection 1. (/"(-^.^et) = 180° ± 15°. UAl algorithm. Li„t = 3 /b"! 



jD+ctust 


30 


20 


15 


10 


5 


ppet 


43.021 


42.771 


42.679 


42.755 


43.202 


p^Jet_ppet 


0.168 


0.167 


0.161 


0.160 


0.127 


Pi Jet 


0.169 


0.168 


0.162 


0.161 


0.128 


^event 


0.033 


0.033 


0.033 


0.033 


0.027 


(m) 


0.100 


0.099 


0.096 


0.099 


0.087 




0.018 


0.018 


0.017 


0.017 


0.014 




4.551 


4.511 


4.470 


4.399 


4.134 




7.054 


6.942 


6.843 


6.777 


6.576 


Nevent(c) 


312191 


287694 


253628 


180811 


40334 


Nevent(b) 


40098 


36223 


30495 


20689 


3639 


29sub/aU 


0.92 


0.91 


0.91 


0.91 


0.90 


Rjet 


0.60 


0.60 


0.60 


0.60 


0.59 


Entries 


56532 


52588 


46991 


34426 


8421 



Table 2: Selection 1. <p(^f.-iet) = 180° ± 15°. LUCELL algorithm. L,„t = 3 fb-^ 



pj-Ctust 
^''CUT 


30 


20 


15 


10 


5 


pfjet 


43.253 


43.000 


42.949 


43.026 


43.408 


p^Jet_ppet 


0.168 


0.165 


0.160 


0.156 


0.121 


ID+Jet 


0.169 


0.166 


0.161 


0.157 


0.121 




0.033 


0.033 


0.033 


0.032 


0.027 




0.103 


0.100 


0.098 


0.093 


0.094 




0.018 


0.018 


0.017 


0.017 


0.015 


p^mtss 


4.556 


4.510 


4.474 


4.382 


4.104 


Tyj-miss 


7.027 


6.915 


6.834 


6.745 


6.595 


Nevent(c) 


304172 


277451 


241228 


164132 


36021 


Nevent(5) 


39256 


34740 


29248 


18937 


3167 


29sub/aU 


0.92 


0.91 


0.91 


0.90 


0.90 


Rjet 


0.65 


0.65 


0.65 


0.64 


0.63 


Entries 


54922 


50723 


44738 


31455 


7751 
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100 < Pt^ < 120 GeV/c 



Table 3: Selection 1. <f>(-yjet) = 180° ± 15°. UAl algorithm. = 3 



ryj-ciust 


30 


20 


15 


10 


5 


Pfjet 


102.627 


104.675 


105.575 


106.329 


106.917 


p^.Jet_ppet 


0.546 


0.538 


0.523 


0.501 


0.488 




0.548 


0.539 


0.525 


0.502 


0.489 


r>v£Jel 
^event 


0.044 


0.042 


0.040 


0.038 


0.034 


Pi Jet 


0.258 


0.249 


0.234 


0.228 


0.216 




0.023 


0.022 


0.021 


0.019 


0.019 


p^miss 


5.166 


5.139 


5.102 


5.053 


4.913 


p^raiss 


14.245 


14.512 


14.817 


14.831 


15.412 


Nevent(c) 


18289 


13417 


10124 


6149 


1051 


Nevent(b) 


2887 


1989 


1435 


779 


113 


29sub/all 


0.89 


0.89 


0.88 


0.88 


0.86 


Rjet 


0.61 


0.61 


0.61 


0.61 


0.60 


Entries 


63316 


48178 


37512 


23472 


4467 



Table 4: Selection 1. <p(.y,jet) = 180° ± 15°. LUCELL algorithm. L^m = 3 /b"^ 



ryxciust 


30 


20 


15 


10 


5 




103.378 


105.266 


106.137 


106.938 


107.216 


Pl"^-Pl"' 


0.549 


0.544 


0.524 


0.475 


0.491 


pfJeL 


0.552 


0.546 


0.525 


0.477 


0.492 


-^event 


0.045 


0.043 


0.041 


0.037 


0.034 




0.260 


0.249 


0.240 


0.223 


0.198 




0.023 


0.022 


0.02 1 


0.019 


0.017 


p-^miss 


5.168 


5.136 


5.110 


5.010 


4.897 


pj.miss 
^^veJet 


14.169 


14.506 


14.527 


14.442 


15.832 


Nevent(c) 


17309 


12498 


9257 


5137 


984 


Neventj-ft-) 


2704 


1866 


1308 


620 


102 


29sub/all 


0.89 


0.89 


0.88 


0.87 


0.85 


Rjet 


0.65 


0.65 


0.65 


0.65 


0.64 


Entries 


59683 


44691 


34139 


20072 


4019 
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200 < Pt^ < 240 GeV/c 



Table 5: Selection 1. <f>(-yjet) = 180° ± 15°. UAl algorithm. = 3 



ryj-ciust 


30 


20 


15 


10 


5 


Pfjet 


211.973 


213.370 


214.124 


214.874 


215.511 


p^.Jet_ppet 


0.886 


0.874 


0.823 


0.768 


0.639 




0.889 


0.877 


0.825 


0.770 


0.640 


r>v£Jel 
^event 


0.039 


0.038 


0.037 


0.035 


0.033 


Pi Jet 


0.420 


0.399 


0.360 


0.338 


0.344 




0.020 


0.019 


0.019 


0.017 


0.016 


p^miss 


5.558 


5.529 


5.444 


5.368 


5.130 


p^raiss 


24.382 


24.517 


23.705 


23.182 


20.955 


Nevent(c) 


1081 


753 


547 


317 


52 


Nevent(b) 


152 


100 


70 


36 


6 


29sub/all 


0.86 


0.85 


0.85 


0.84 


0.82 


Rjet 


0.62 


0.62 


0.62 


0.62 


0.61 


Entries 


52542 


37741 


28477 


17189 


3142 



Table 6: Selection 1. 4>(-y,jet) = 180° ± 15°. LUCELL algorithm. L^nt = 3 fb~^ 



ryxciust 


30 


20 


15 


10 


5 




212.521 


213.982 


214.736 


215.460 


216.044 


Pl"^-Pl"' 


0.866 


0.850 


0.802 


0.742 


0.568 


pfJeL 


0.869 


0.853 


0.805 


0.744 


0.569 


-^event 


0.038 


0.037 


0.036 


0.034 


0.028 




0.417 


0.390 


0.353 


0.336 


0.268 




0.019 


0.019 


0.019 


0.017 


0.014 


p-^miss 


5.529 


5.487 


5.412 


5.337 


4.975 


pj.miss 
^^veJet 


24.076 


24.016 


23.622 


23.102 


21.347 


Nevent(c) 


1012 


694 


487 


261 


44 


Neventj-ft-) 


138 


90 


60 


27 


4 


29sub/all 


0.86 


0.85 


0.84 


0.83 


0.80 


Rjet 


0.66 


0.65 


0.65 


0.65 


0.64 


Entries 


49253 


34775 


25582 


14562 


2786 
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300 < Pt^ < 360 GeV/c 



Table?: Selection 1. <f>(-yjet) = 180° ± 15°. UAl algorithm. = 3 



ryj-ciust 


30 


20 


15 


10 


5 


Pfjet 


320.158 


321.502 


322.289 


322.869 


322.911 


p^.Jet_ppet 


1.077 


1.069 


1.060 


1.015 


1.161 




1.081 


1.072 


1.063 


1.018 


1.163 


r>v£Jel 
^event 


0.035 


0.034 


0.033 


0.033 


0.035 


Pi Jet 


0.515 


0.506 


0.476 


0.448 


0.433 




0.019 


0.018 


0.017 


0.017 


0.017 


p^miss 


5.764 


5.721 


5.692 


5.572 


5.691 


p^raiss 


31.983 


32.597 


33.078 


32.371 


35.201 


Nevent(c) 


172 


117 


84 


46 


8 


Nevent(b) 


25 


15 


10 


6 


1 


29sub/all 


0.84 


0.83 


0.82 


0.80 


0.78 


Rjet 


0.62 


0.62 


0.62 


0.62 


0.61 


Entries 


46297 


32513 


24157 


14318 


2642 



Table 8: Selection 1. 4>(-y,jet) = 180° ± 15°. LUCELL algorithm. L^nt = 3 fb~^ 



ryxciust 


30 


20 


15 


10 


5 




320.687 


322.011 


322.732 


323.248 


323.646 


Pl"^-Pl"' 


1 .072 


1.061 


1.055 


0.983 


1.052 


pfJeL 


1.076 


1.064 


1.057 


0.985 


1.053 


-^event 


0.035 


0.033 


0.033 


0.032 


0.032 




0.507 


0.483 


0.480 


0.412 


0.388 




0.018 


0.018 


0.017 


0.016 


0.013 


p-^miss 


5.761 


5.722 


5.686 


5.499 


5.465 


pj.miss 
^^veJet 


32.220 


33.054 


33.221 


31.968 


34.611 


Nevent(c) 


161 


106 


74 


39 


7 


Neventj-ft-) 


22 


14 


9 


5 


1 


29sub/all 


0.83 


0.82 


0.81 


0.80 


0.77 


Rjet 


0.66 


0.65 


0.65 


0.65 


0.64 


Entries 


43320 


29783 


21707 


12104 


2334 
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